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ABSTRACT 

We studied the outburst properties of the hyper-luminous X-ray source ESO 243-49 HLX-1, using the 
full set of Swift monitoring observations. We quantified the increase in the waiting time, recurrence 
time, and e-folding rise timescale along the outburst sequence, and the corresponding decrease in out¬ 
burst duration, total radiated energy, and e-folding decay timescale, which confirms previous findings. 

HLX-1 spends less and less time in outburst and more and more time in quiescence, but its peak 
luminosity remains approximately constant. We compared the HLX-1 outburst properties with those 
of bright Galactic low-mass X-ray binary transients (LMXBTs). Our spectral analysis strengthens 
the similarity between state transitions in HLX-1 and those in Galactic LMXBTs. We also found that 
HLX-1 follows the nearly linear correlations between the hard-to-soft state transition luminosity and 
the peak luminosity, and between the rate of change of X-ray luminosity during the rise phase and 
the peak luminosity, which indicates that the occurrence of the hard-to-soft state transition of HLX-1 
is similar to those of Galactic LMXBTs during outbursts. We found that HLX-1 does not follow the 
correlations between total radiated energy and peak luminosity, and between total radiated energy 
and e-folding rise/decay timescales we had previously identified in Galactic LMXBTs. HLX-1 would 
follow those correlations if the distance were several hundreds of kiloparsecs. However, invoking a 
much closer distance for HLX-1 is not a viable solution to this problem, as it introduces other, more 
serious inconsistencies with the observations. 

Keywords: accretion, accretion disks — black hole physics —X-rays:binaries — X-rays: individual 
(HLX-1) 


1. INTRODUCTION 

HLX-1 is the most luminous ultraluminous X-ray 
source (ULX; Feng & Soria 2011, for a review) known 
with a maximum isotropic X-ray luminosity of > 
10"^^ erg s“^ (Farrell et al. 2009; Godet et al. 2009). 
The extremely high luminosity makes HLX-1 currently 
the best candidate to harbor an intermediate-mass black 
hole (IMBH, 10^-10® Mq) among ULXs. Although other 
ULXs in nearby galaxies were initially flagged as can¬ 
didate IMBHs, recent works show that some of them 
contain a stellar-mass black hole (BH; Liu et al. 2013; 
Motch et al. 2014) or even a neutron star (NS; Bachetti 
et al. 2014). HLX-1 is projected in the sky within the 
D25 (the diameter of the light isophote having a surface 
brightness of 25.0 mag arcsec“^) of the SOa spiral galaxy 
ESO 243-49 (redshift distance «95 Mpc), « 8" north¬ 
east of its nucleus (Farrell et al. 2009). The redshift 
measurement of the Ha emission line from the optical 
counterpart of HLX-1 confirmed the association of HLX- 
1 with ESO 243-49 (Wiersema et al. 2010; Soria et al. 
2013). However, there is an offset of «430 km s“^ be¬ 
tween the recession velocity of HLX-1 and that of the 
nucleus of ESO 243-49: this suggests that HLX-1 may be 
in a dwarf satellite galaxy or star cluster near ESO 243-49 
rather than in the galaxy itself (Soria et al. 2013). 
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Monitoring observations of HLX-1 with Swift/XKT 
have shown that the source undergoes recurrent out¬ 
bursts (almost once a year), where the outburst pro¬ 
files are well described by a fast-rise-exponential-decay 
(FRED) form (see also Lasota et al. 2011; Godet et al. 
2014), similar to the outbursts of many bright Galactic 
low-mass X-ray binary transients (LMXBTs; see Chen 
et al. 1997; Yan & Yu 2015). During its outbursts, HLX- 
1 displayed spectral state transitions which also resemble 
those of the bright Galactic LMXBTs (Godet et al. 2009; 
Servillat et al. 2011). Webb et al. (2012) have detected 
transient radio emission following the hard-to-soft state 
transition during the 2010 and 2011 outbursts. By anal¬ 
ogy with the radio flares seen in Galactic BH LMXBs 
(Fender et al. 2004), the transient radio emission was in¬ 
terpreted as discrete jet ejection events during the state 
transitions (Webb et al. 2012). 

It is generally believed that the outbursts of bright 
Galactic LMXBTs are triggered by the thermal-viscous 
disk instability (Lasota 2001, for a review). In this 
model, the outburst duration is approximately the vis¬ 
cous timescale at the minimum unstable disk radius, 
which is the radius at which the irradiated disk becomes 
cold enough to be dominated by neutral hydrogen. If 
HLX-I is located at «95 Mpc, then its maximum lumi¬ 
nosity, « 10'^^ erg s“^, implies a viscous timescale ~100 
years (see Lasota et al. 2011, 2015, assuming that the 
disk is large enough to allow for a neutral outer region). 
However, the observed outburst duration of HLX-1 is 
^100 days, i.e., at least 2 orders of magnitude shorter 
than needed, therefore suggesting that the outbursts are 
triggered from a much smaller radius than required by 
the thermal-viscous disk instability model. Clearly, an- 
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other mechanism is probably needed to trigger the out¬ 
bursts. For example, Lasota et al. (2011) proposed that 
outbursts are triggered by the passage at periapse of the 
companion star on an eccentric orbit. However, this 
model cannot easily explain why the outburst repeats 
at irregular intervals (between «370 and 460 days: Kong 
2011; Lasota et al. 2011; Godet et al. 2013, 2014, 2015; 
Miller et al. 2014; Kong et al. 2015). Based on the ob¬ 
served decline of the peak count rate, integrated lumi¬ 
nosity, duration, and decay time over the sequence of 
outbursts, it was also suggested that the donor star may 
be undergoing tidal disruption (Godet et al. 2013, 2014; 
Miller et al. 2014). Miller et al. (2014) proposed an al¬ 
ternative model based on (low-angular-momentum) wind 
accretion instead of Roche lobe overflow. King & Lasota 
(2014) suggested that HLX-1 may be a stellar-mass bi¬ 
nary system like SS433 (see also Lasota et al. 2015), 
in which the X-ray emission comes from the beamed jet. 
None of those scenarios appear to be fully consistent with 
all the observational constraints, and there is still no gen¬ 
eral consensus on what causes the outbursts. 

In previous work (Yan & Yu 2015), we have system¬ 
atically studied the outburst properties of bright Galac¬ 
tic LMXBTs, quantifying statistical correlations between 
observable quantities and comparing them with the pre¬ 
dictions of the thermal-viscous disk instability model. 
In order to understand the nature of HLX-I outbursts, 
in this paper, we use the same parameters as used by 
Yan & Yu (2015) to quantify the outburst properties of 
HLX-I, how they changed from outburst to outburst, 
and how they compare with the observed properties of 
bright Galactic LMXBTs outbursts. 

2. DATA ANALYSIS AND RESULTS 
2.1. Different Spectral States 

The Swift X-Ray Telescope (XRT) has monitored 
HLX-I regularly for more than six years, providing us 
with an opportunity to investigate the X-ray temporal 
and spectral variability of this source. We obtained the 
0.3-10 keV count rate and the (1.5-10 keV)/(0.3-1.5 
keV) hardness ratio (HR) of each Swift/XRT observation 
from the online Swift/XKT product generator (Evans 
et al. 2007, 2009). We plotted the HR distribution from 
all observations with positive HR (Figure 1): it shows 
two peaks, suggesting that HLX-1 has two distinct spec¬ 
tral states (hard state and soft state) by analogy with 
Galactic X-ray binaries (XRBs; Yu & Yan 2009; Yan & 
Yu 2010). The presence of different spectral states in 
HLX-1 was first noted by Godet et al. (2009) and was 
also investigated by Servillat et al. (2011). We used the 
double-peaked HR distribution to define the hard and 
soft states: we htted the histogram with two Gaussian 
components (Figure 1) and found that the centroids are 
roughly at HR « 0.1 and HR « 0.6. We classify all the 
observations with HR < 0.1 as in the soft state, and all 
those with HR > 0.6 as in the hard state. The obser¬ 
vations corresponding to 0.1 < HR < 0.6 consist of a 
mixture of soft, hard and perhaps intermediate states: 
we classify those observations as the transitional state. 

Hardness—intensity diagrams (HIDs) are a very useful 
tool to characterize the different spectral states in Galac¬ 
tic LMXBs {e.g. Miyamoto et al. 1995; Homan et al. 
2001; Maccarone & Coppi 2003; Kording et al. 2008; 



Figure 1. Histogram of the (1.5-10 keV)/(0.3-1.5 keV) HR dis¬ 
tribution (solid line), fit with a double Gaussian model (dotted 
line). The Gaussian peaks are roughly at HR Ri 0.1 and HR ^ 0.6. 
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Figure 2. Srof/i/XRT hardness—intensity diagram of HLX-1. 
We include all observations with detection significance > 2cr in 
the full 0.3—10 keV band (Evans et al. 2007, 2009). For clarity, 
we only plot two points with error bars to indicate representative 
uncertainty levels. The dotted lines separate the soft, transitional 
and hard spectral states. For each outburst, dashed lines indicate 
the spectral evolution from the hard to the soft state and then back 
to the hard state. 


Belloni 2010). In Figure 2 we plot the HID of all HLX- 
I outbursts. The diagram clearly shows two clusters of 
data points, corresponding to the high/soft and low/hard 
state. In each outburst, HLX-1 evolved from the hard 
state to the soft state during the rise phase and then re¬ 
turned to the hard state during the decay phase. This be¬ 
havior is similar to most outbursts of Galactic LMXBTs 
{e.g. Homan et al. 2001; Belloni et al. 2006; Belloni 2010; 
Dunn et al. 2010; Munoz-Darias et al. 2014). 

We extracted stacked spectra of the Swift/XKT obser¬ 
vations in each of the three states (Figure 3), using the 
online XRT product generator (Evans et al. 2009), which 
also provides the associated background and response 
files. We only included observations with detection sig¬ 
nificance > 2(7. We analyzed and modeled the spectra 
with XSPEC Version 12.8.2 (Arnaud 1996). Because of 
the low count rate, we grouped the spectrum of the hard 
state to at least 1 count bin“^, and applied the Gash 
statistics (Cash 1979), and grouped the spectra of the 
transitional and soft sates to at least 20 counts bin“^, so 
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Table 1 

Best-fitting parameters in different spectral states 


Spectral States 

Tin (keV) 

-^disk 

r 

ATpow (10-5) 

X^/dof 

Hard State 

Transitional State 

Soft State 

0 16+0-02 

0 

^•^^-0.01 

25.60t3|83 

32.08t7;i2 

-L-0^_0.32 

2 45+0.41 

^•^^-0.45 

^•OO_0_g2 

1.DU_q 3g 

1 62+^’^^ 

9 1 q+2.20 

^•^'^-1.33 

67.27/104 (C-statistic) 
38.26/52 

81.52/82 



Figure 3. Stacked Swift/X.KT spectra and best-fitting models for 
the three states defined in Section 2.1: red data points for the soft 
state, green for the transitional state, and blue for the hard state. 
See Table 1 for the best-fitting parameters. 


that we could use statistics. We ignored photon ener¬ 
gies below 0.5 keV. We have tried to free the column den¬ 
sity, but it can only be well constrained by the spectrum 
of the soft sate at a value of A^h = x 10^° cm“^, which 
is consistent with the value used before (Farrell et al. 
2009; Webb et al. 2010, 2012). Therefore, we fixed the 
column density at the line-of-sight value TVh = 5 x 10^° 
cm“^ for all the three spectra. 

The stacked spectrum of the hard state is well fit by a 
single power-law model with photon index T « 1.6 (Ta¬ 
ble 1), similar to the typical values found in the hard 
state of LMXBs. The stacked spectrum in the soft state 
is fit by a disk-blackbody model plus a power-law compo¬ 
nent that accounts for significant residual emission above 
2 keV. The disk contributes «85% of the 0.3-10 keV flux, 
similar to the situation in the disk-dominated soft state of 
BH LMXBs. The peak color temperature Tin in the soft 
state is r:! 0.2 keV (Table 1), which is much less than the 
typical value in luminous BH LMXBs (Tin ~ 1 keV). The 
corresponding inner-disk radius i?in 50, 000- 100, 000 
km for a plausible range (10°-80°) of disk inclination an¬ 
gles, at a distance of 95 Mpc. The low color temperature 
and large radius of the thermal emitter are consistent 
with a disk around an IMBH (Godet et al. 2009; Davis 
et al. 2011; Servillat et al. 2011). However, the ther¬ 
mal X-ray emission components with cool temperatures 
and large radii are also seen in ultraluminous supersoft 
sources where the emission is probably coming from the 
photosphere of a thick outflow rather than from a disk 
(Poutanen et al. 2007; Shen et al. 2015), and therefore 


do not require an IMBH as the accreting object {e.g. the 
supersoft ULX in M 101 is probably powered by a stellar- 
mass BH, Liu et al. 2013). Finally, the stacked spectra 
of HLX-1 in the transitional state are well fit by a disk- 
blackbody (Tin ~ 0.16 keV) plus power-law (F « 2.45) 
model (Table 1), where the disk emission accounts for 
«52% of the 0.3-10 keV flux. The disk normalization 
is similar to the soft state (Table 1), which is consistent 
with a disk extending down to the inner-most stable cir¬ 
cular orbit (ISCO). 

2.2. State Transition Luminosity 

According to our definition of the three spectral states, 
the mean XRT count rates are «0.022 ct s“^, «0.005 ct 
s“^ and «0.012 ct s“^ for the soft, hard and transitional 
states, respectively. Based on our spectral fitting results 
(Table 1), the unabsorbed 0.3-10 keV fluxes correspond¬ 
ing to an XRT count rate of 1 ct s“^ are 3.11 x 10“^^ erg 
s“^ cm“^ in the soft state, 2.60 x 10“^^ erg s“^ cm“^ in 
the hard state, and 3.03 x 10“^^ erg s“^ cm“^ in the tran¬ 
sitional state. For each observation, we used the hard¬ 
ness ratio to identify the spectral state, then converted 
the count rate to an unabsorbed flux by using the appro¬ 
priate conversion factor. The ^wi/t/XRT light curve in 
units of erg cm is shown in Figure 4. 

The Swift observations have so far covered seven out¬ 
bursts of HLX-1, plus only a few data points for an out¬ 
burst in 2008 (Figure 4); we exclude those 2008 data 
points in the rest of this paper. We determined the 
luminosity Thard-to-soft at the transition from hard to 
soft state during the outburst rise, and the luminosity 
Tsoft-to-hard at the reverse transition during the decay 
phase, using the same method previously applied by Yu 
& Yan (2009) to Galactic LMXBTs. In that method, we 
used the luminosity of the last hard state data point as 
the hard-to-soft state transition luminosity (Yu & Yan 
2009; Tang et al. 2011). However, it is much harder to 
catch the actual transition in the HLX-1 outbursts be¬ 
cause the rise time is short compared with the Swift sam¬ 
pling frequency. Therefore, we note that we might be un¬ 
derestimating Thard-to-soft with our method. The Swift 
observations do not cover the hard-to-soft state transi¬ 
tion in the 2009 and 2010 outbursts; for the other four 
outbursts, we list Thard-to-soft in Table 2. The first ob¬ 
servation with HR > 0.6 during the decay phase was 
identified as the first hard state observation; we used the 
luminosity of this observation as the soft-to-hard state 
transition luminosity Tgoft-to-hard (Yu & Yan 2009). We 
list the Tsoft-to-hard for six outbursts in Table 2. 

We had already found a correlation between the hard- 
to-soft state transition luminosity and the peak lumi¬ 
nosity in Galactic XRBs (Yu et al. 2004; Yu & Dolence 
2007; Yu & Yan 2009; Tang et al. 2011). We fit the cor- 
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Figure 4. Complete Swift/XRT light curve (unabsorbed flux in unit of erg s ^ cm ^), binned to individual observations. 



Figure 5. Correlation between ihard-to—soft ^peak in Galac¬ 
tic LMXBTs, and the location of HLX-1 outbursts (plotted as tri¬ 
angles) in the same parameter space. The solid line is a linear 
regression of the form logLpeak = A -|- B X logLhard-to-soft, ht 
to the Galactic LMXBT sample using the method of Kelly (2007). 
The dashed lines show the intrinsic scatter and the shaded region 
shows the Icr confidence interval of the regression line. HLX-1 is 
roughly consistent with the correlation. 


relation in Galactic XRBs with a function log Lpeak = 
A + B X logLhard-to-soft by using the Bayesian approach 
in Kelly (2007), and obtained a best-fitting slope of B 
= 0.96±0.04, which means that Lhard-to-soft is linearly 
proportional to Tpeak- Now we have added the values 
of Lhard-to-soft and Lpeak for the four HLX-1 outbursts 
(Figure 5). HLX-1 roughly follows this correlation we 
found in Galactic XRBs, which confirms the prediction 
that the correlation extends to ULXs (Yu & Yan 2009). 

2.3. Outburst Parameters 


Table 2 

State Transition Luminosity 


Outburst 

-^hard —to —soft 

(10^^ erg s“^) 

-^soft —to —hard 

(10"^^ erg s“^) 

2009 


0.62 ±0.44 

2010 


1.25 ±0.29 

2011 

0.90 ±0.51 

2.32 ± 1.35 

2012 

2.74 ± 1.32 

0.55 ±0.48 

2013 

2.30 ± 1.04 

0.90 ±0.44 

2015 

1.10 ±0.62 

1.78 ±0.69 


The shape of the X-ray light curves for the six HLX- 
1 outbursts with good Swift coverage are qualitatively 
similar. The 2009 and 2011 outbursts show the clearest 
FRED profiles (Figure 6), which is typical of Galactic 
LMXBTs outbursts (Chen et al. 1997). We only have 
flux upper limits for the early rise phase of the 2010 out¬ 
burst, but they are also consistent with a FRED profile 
(Soria 2013, Fig. 1). The 2012, 2013, and 2015 outbursts 
are roughly consistent with FRED profiles, but also show 
secondary peaks or reflarings during the decline (Fig¬ 
ure 6). 

We measured the main outburst parameters for the 
six outbursts (Table 3, Figure 7), including peak lumi¬ 
nosity, e-folding rise and decay timescales, rate of change 
of the luminosity during the rise and decay, outburst du¬ 
ration, and total radiated energy, according to the defi¬ 
nitions in Yan & Yu (2015). The peak luminosity Lpeak 
is defined as the maximum luminosity in each outburst. 
We defined L,.ise^io%, T(iecay, 10 %; -brise,90% ffod L(iecay,90% 
as the measured data points closest to 10% and 90% 
of Lpeak during the rise and decay phases, respectively. 
The e-folding rise timescale is defined as rnse = At/lnT, 
where At is the time interval between the data points 
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Figure 6. Swift/XKH' light curves of each outburst (excluding the few sparse data points in 2008). In each panel, Time = 0 corresponds 
to the time of Trise io%- 


Table 3 

Outburst parameters 


Peak Time 

(MJD) 

T^peak 

(xl0‘‘2) 

(erg s-i) 

Duration 

(days) 

'-p a 

waiting 

(days) 

T1 ^ 

J- J-recurrence 

(days) 

T9 ^ 

■^recurrence 

(days) 

E 

(xlO"^®) 

(erg) 

"^rise 

(days) 

"^decay 

(days) 

-^rise 

(xlO^^) 

(erg day“^) 

-^decay 

(xlO^O) 

(erg day“^) 

55060 

1.2±0.1 

188 

<265 

<454 

>286 

9.7±0.1 

3±0 

62±16 

11.3±0.9 

0.6±0.1 

55446 

1.3±0.2 

128±8 

192±8 

320 

380 

6.6±0.4 


48±6 


1.0±0.2 

55792 

1.4±0.4 

110 

224 

334 

352 

5.9±0.2 

2±1 

50±17 

17.2±5.9 

1.1±0.3 

56178 

1.2±0.2 

96 

264 

360 

375 

4.8±0.2 

6±2 

25±8 

7.0±1.5 

1.4±0.3 

56586 

l.OiO.l 

84 

308 

392 

404 

4.3±0.1 

11±3 

18±4 

3.4±0.6 

2.0±0.4 

57060 

1.4±0.8 

69 

386 

455 

470 

3.7±0.1 

14±5 

18±6 

4.1±2.8 

3.0±2.1 


^ Waiting time 


Recurrence time defined as the interval between .^{jgcay 10% outburst and 10 % the previous outburst 

Recurrence time defined as the interval between i'pjgg 10% outburst and i'j.jgg io% the previous outburst 

Total radiated energy per outburst 


corresponding to ii.ise,io% and irise, 90 %! and A is the ra¬ 
tio of l^rise, 90 %/-^rise,io%- T^e e-folding decay timescale 
Tdecay is defined in the same way, but for i(jecay,io% and 

idecay, 90 %- The rate of change of luminosity Lrise and 
idecay are defined as (irise.90%-.^rise.io%)/^t for the rise 
phase and as (Ldecay, 90 % - idecay.io%)/At for the decay 
phase, respectively. The outburst duration is defined as 
the time interval between Trise,io% and T(jecay,io%- We 
have two definitions of the recurrence time, one is the in¬ 
terval between idecay,io% of an outburst and Tdecay,io% 
of the previous outburst (abbreviated as T1 recurrence); 
the other is the interval between Trise,io% of an outburst 
and Trise,io% of the previous outburst (abbreviated as 
T2i.ecurrence). The total radiated energy E is defined as 
the integrated luminosity from the time of Trise,io% to 
the time of Tdecay,io% (for this integration, the data gaps 


due to sparse coverage were filled with linear interpola¬ 
tions). Finally, we determined the waiting time for each 
outburst, defined as the time interval between Trise,io% 
of an outburst and Tdecay,io% of the previous outburst. 
For the outburst in 2009, we only have upper limits to 
the waiting time and Tlrecurrence and lower limits to the 
T2recurrence, because the data coverage of the 2008 out¬ 
burst is too limited to permit the identification of the 
times of Lj-ise^ioyo and Tdecay,10%- 

The uncertainties in the outburst durations, wait¬ 
ing times, and recurrence times are mainly due to the 
Swift/XBT observation sampling frequency. The typical 
interval of observations in different outbursts is in the 
range of 1-14 days. For example, for the 2010 outburst 
we are not able to identify Trise, 90 % (Figure 6); therefore, 
we cannot calculate the Trige and Lrise- Instead, we use 
the earliest significant detection (Figure 6) as Trise,io% 
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Outburst squence 

Figure 7. Evolution of the main outburst parameters defined in 
Section 2.3. The total radiated energy E, the outburst duration, 
and the e-folding decay timescale rdecay decreased and the waiting 
time, recurrence time, and e-folding rise timescale Tj-[se increased. 
The first outburst is the one in 2009. Twaiting is the waiting time. 
Tlrecurrence is the recurrence time defined as the interval between 
Ldecay,io% outburst and Tdecay,io% previous outburst. 

T2recurrence is the recurrence time defined as the interval between 
Lrise 10% outburst and of the previous outburst 


to estimate lower limits to the duration, total radiated 
energy, and waiting time of the 2010 outburst. The ob¬ 
servation before the earliest significant detection is only 
an upper limit: by taking that value as Trise,io%i esti¬ 
mated upper limits to the duration, total radiated energy, 
and waiting time (Table 3). 

Our main hndings are the following. 

1. There is no clear trend for the peak luminosity 
with outburst sequence, probably because it is well 
constrained in the range of 1.0-1.4 xlO"'^ erg s“^ 
(Figure 7, Table 3). The value and uncertainty of 
the peak luminosity of the 2015 outburst is quite 


large, which could be due to the short-term vari¬ 
ability, since the exposure time of this observation 
is short. Godet et al. (2013, 2014) argued that the 
peak count rate decreases with outburst sequence, 
but we did not find such a trend in the peak lu¬ 
minosity. There may be two reasons for this dis¬ 
crepant conclusion. Firstly, our light curve was re¬ 
binned by observation, while Godet et al. (2013, 

2014) rebinned theirs by photon count (20 counts 
per bin): as a result, the peak values are slightly 
different. Secondly, our light curve is in flux rather 
than count rate units: the conversion factor is dif¬ 
ferent depending on whether the peak of the out¬ 
burst was reached during the transitional or soft 
state. 

2. The e-folding rise timescale increases at least since 
the 2011 outburst. 

3. The e-folding decay timescale decreases, which con¬ 
firms the result in Miller et al. (2014). 

4. The outburst duration also decreases along the out¬ 
burst sequence , which confirms previous results in 
Godet et al. (2013, 2014), from «188 days in the 
2009 outburst to «69 days in the 2015 outburst 
(Figure 7, Table 3). 

5. The total radiated energy decreases along the out¬ 
burst sequence, which confirms the result in Miller 
et al. (2014), from E = 9.7 ± 0.1 x 10^® erg in the 
2009 outburst to E = 3.7 ± 0.1 x 10"*® erg in the 
2015 outburst (Figure 7 and Table 3). This corre¬ 
sponds to a mass transfer of «5 x 1 O“®M 0 in the 
2009 outburst, down to «2 x 10~^Mq in the 2015 
outburst (assuming a standard efhciency wO.l). 

6 . The waiting time increases along the outburst se¬ 
quence at least since the 2010 outburst, from «192 
days for the 2010 outburst to wdSb days for the 
2015 outburst (Figure 7, Table 3). 

7. The recurrence time increases along the outburst 
sequence, from « 320 days of the 2010 outburst to 
« 455 days of the 2015 outburst for Tlrecurrence, 
and from « 352 days of the 2011 outburst to « 
470 days of the 2015 outburst for T2recurrence (Fig¬ 
ure 7, Table 3) . The increase of the recurrence time 
creates serious problems for its original and sim¬ 
plest interpretation as the binary period. There¬ 
fore, Godet et al. (2014) tried to resolve this prob¬ 
lem by considering a binary system with extreme 
parameters, and they predicted an increase of the 
binary period in such system, which is consistent 
with the increase of recurrence time (Godet et al. 

2015) . 

8 . The duty cycle (defined as the fraction of time 
spent in outbursts) is «0.3. 

2.4. Comparison of Outburst Properties with Bright 
Galactic LMXBTs 

The comparison between the outburst properties of 
HLX-1 and those we determined previously for bright 
Galactic LMXBTs (Yan & Yu 2015) can be summarized 
as follows. 
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Figure 8. Relations between L and Lpeak during outburst rise 
(top panel) and decay (bottom panel) in Galactic LMXBTs and 
HLX-1. The solid line shows the best-fitting linear correlation for 
the Galactic sources, with their intrinsic scatter (dashed lines) and 
l(j confidence intervals (shaded regions). HLX-1 follows the cor¬ 
relation between Lrise and Lpeak? but not the one between Ldecay 
and Lpeak- HLX-1 would have to be located at p=; 540 kpc to satisfy 
the both correlations, which is inconsistent with its optical redshift 
(see details in Section 3.3). 



Figure 9. Relation between E and Lpeak Galactic LMXBTs 
and HLX-1. The solid line shows the oest-fitting correlation for 
the Galactic sources, with their intrinsic scatter (dashed lines) and 
Icr confidence intervals (shaded regions). HLX-1 does not follow 
this correlation. As for the other cases shown in Figure 8, HLX-1 
will follow this correlation if it is located at a distance of 5^540 kpc 
(see details in Section 3.3), which is inconsistent with its optical 
redshift. 



Figure 10. Relation between e-folding rise (top panel) and de¬ 
cay (bottom panel) timescales r and total radiated energy E, for 
Galactic LMXBTs and HLX-1. Solid and dashed lines, and shaded 
regions, are defined as in Figures 8 and 9. HLX-1 is several orders 
of magnitude off the correlation. In this case, too, it would sat¬ 
isfy the LMXBT correlation if it were at 5^540 kpc (see details in 
Section 3.3). 


1. The average peak luminosity of HLX-1 (« 1.2 x 
10^^ erg s“^) is more than four orders of magni¬ 
tude larger than the average peak luminosity of 
bright Galactic LMXBTs ( « 7 x 10^^ erg s“^ and 
« 4 X 10^^ erg s“^ for BH and NS LMXBTs, re¬ 
spectively) . 

2. The average total radiated energy per outburst of 
HLX-1 (« 5.8 X 10^® erg) is about five orders of 
magnitude larger than those measured in bright 
Galactic LMXBTs ( « 2 x 10'^'^ erg and « 5 x 10^^ 
erg for BH and NS LMXBTs, respectively). 

3. The average outburst duration of HLX-1 (« 113 
days) is about a factor of 3 longer than that of the 
NS LMXBTs outbursts (« 39 days), and is a little 
bit longer than the average outburst duration of 
BH LMXBTs (« 88 days), but still within the Icr 
range of BH LMXBTs. 

4. The average e-folding rise timescale of HLX-1 («7 
days) is about two times larger than the Galac¬ 
tic NS LMXBTs («4 days), and is similar to the 
Galactic BH LMXBTs («8 days). 

5. The average e-folding decay timescale of HLX-1 
(r:! 37 days) is a few times longer than for Galac¬ 
tic NS LMXBTs (average of «!! days) and BH 














































































































LMXBTs (average of !^26 days), but still within 
the Icr range of values measured for BH LMXBTs. 

Positive correlations have been found between the rate 
of change of luminosity and the peak luminosity in both 
the rise and decay phases, for Galactic LMXBTs (Yan 
& Yu 2015). Using a Bayesian linear regression (Kelly 

2007), we had found that loginse = A + B x logLpeak 
with B = 0.99±0.20, and logLdecay = + X log Lpeak 

with B' = 0.62 ± 0.16. We extrapolated those empirical 
relations to see whether HLX-1 is also consistent with 
them (Figure 8). We found that HLX-1 follows the corre¬ 
lation between Lrise and Lpeak, but not the one between 
idecay and Lpeak (Lc., its decay timescale is too short 
to be a scaled-up LMXBT; see also Lasota et al. 2015). 
Other correlations found in Galactic LMXBTs are those 
between the total radiated energy and peak luminosity, 
and those between the e-folding rise/decay timescales 
and total radiated energy; none of those relations is sat¬ 
isfied by HLX-1 (Figure 9 and Figure 10). (Probably 
only by coincidence, HLX-1 would satisfy those Galactic 
LMXBTs correlations if it were located at a distance of 
^540 kpc, which is inconsistent with its optical line red- 
shift and is implausible based on other arguments, see 
Section 3.3.). 

2.5. Comparison with BH LMXB T HI 74 3- 322 

In addition to comparing the outburst properties of 
HLX-1 with the average outburst properties of bright 
Galactic LMXBTs (Section 2.4), we compare HLX-1 
with some individual systems which undergo regular out¬ 
bursts. For example, the BH LMXBT H1743—322 has 
been in outburst nearly once a year for the last 10 years 
(Figure 11), as shown by the Swift/BP^A light curve 
(Krimm et al. 2013). We used that light curve to quantify 
the outburst properties of H1743—322, and investigate 
its similarities with HLX-1. We only considered those 
outbursts with a peak count rate larger than 0.1 crab 
at a 2(7 significance level; this means all outbursts after 
MJD 54700 (j.e., after 2008 August). We used the same 
method described in Section 2.3 to measure the outburst 
parameters. The outbursts occur quite regularly with an 
average recurrence time of 274 ± 56 days, which is very 
similar to that of HLX-1 (Figure 4). We found that the 
waiting time and recurrence time first decrease and then 
increase (see Figure 12), possibly consistent with the be¬ 
havior of HLX-1 (bottom panels of Figure 7). However, 
the waiting time and recurrence time suddenly decrease 
to the average level in the latest 2015 outburst Figure 12. 
Instead, we did not find any clear trend in the outburst 
duration, fluence, peak flux, and e-folding rise and decay 
timescales versus outburst sequence. 

We repeated the outburst analysis for H1743—322 us¬ 
ing RXTE/ASM- and MAXI data (2-10 keV; Matsuoka 
et al. 2009). We obtained similar values for the wait¬ 
ing and recurrence times; this shows that such parame¬ 
ters are not significantly affected by the choice of energy 
band. We do not have enough evidence at this stage to 
test whether or not the possible similarity in the out¬ 
burst waiting and recurrence times between HLX-1 and 
H1743—322 points to the same physical mechanism for 
outburst triggering and evolution. 

For completeness, we also performed a similar compar¬ 


ison between the BH LMXBT GX 339-4, 4U 1630-472, 
the NS LMXBT Aql X-1 (using i?ArA/ASM-pMAXI 
data), and HLX-1. These three LMXBTs also show 
frequent outbursts. However, for GX 339—4 and 4U 
1630—472, we do not find any clear trends when using 
data from their past 5 outbursts and 10 outbursts, nor 
we find that for Aql X-1 when using data from its past 
14 outbursts. 

3. DISCUSSION 

3.1. BH Mass Estimation Based on Lsoft-to-hard 

In Galactic XRBs, the transition luminosity from the 
soft to the hard state during the decay phase spans a 
very narrow range of values ~ l%-4% LEdd (Maccarone 
2003; Kalemci et al. 2013). The average soft-to-hard 
state transition luminosities of HLX-1 is « 1.2 x 10"^^ 
erg s“^ in the 0.3-10 keV band, with a standard devia¬ 
tion of 0.6 X 10"*^ erg s“^. We used the same bolometric 
correction in Maccarone (2003) to convert from 0.3-10 
keV to bolometric luminosity. If the soft-to-hard state 
transition luminosity of HLX-1 is at 2% ^Edd by analogy 
to the Galactic XRBs, the mass of the accreting compact 
object is (8 ± 4) x This is consistent with the 

IMBH mass of ~1O'^-1O®M0 derived from detailed X-ray 
spectral modelling (Farrell et al. 2010; Davis et al. 2011; 
Servillat et al. 2011; Godet et al. 2012; Webb et al. 2012). 

3.2. Sub-Eddington State Transitions or 
Super-Eddington Geometric Beaming? 

The crucial unsolved problem is whether the spectral 
state evolution supports the sub-Eddington IMBH sce¬ 
nario or the super-Eddington stellar-mass BH scenario. 
We address this question in two steps: firstly, we test 
whether our results are consistent with an IMBH, sec¬ 
ond, we discuss whether the super-Eddington scenario 
provides a more plausible explanation. 

Our spectral analysis confirmed that HLX-1 exhibits 
two distinct X-ray spectral states: the low/hard and 
high/soft states (Godet et al. 2009; Servillat et al. 2011). 
In each outburst, HLX-1 evolved from the hard state to 
the soft state during the rise phase and then returned 
to the hard state during the decay phase (Figure 2). 
This spectral state evolution is similar to that of Galac¬ 
tic XRBs. HLX-1 follows the correlation between the 
hard-to-soft transition luminosity and the outburst peak 
luminosity which we have found in Galactic XRBs (Fig¬ 
ure 5). This correlation, which spans six orders of magni¬ 
tude, strongly suggests that the same mechanism deter¬ 
mines the hard-to-soft state transition luminosity during 
an outburst for both Galactic XRBs and HLX-1. The 
X-ray spectrum in the hard state of HLX-1 is dominated 
by a power-law component with a photon index of «1.6 
(Table 1), which is similar to that of Galactic XRBs in 
the hard state. The X-ray spectrum in the soft state is 
dominated by a thermal component (Table 1), which is 
also similar to that of the soft state of Galactic XRBs but 
with lower temperature and larger radius, as expected for 
an IMBH. The luminosity and peak color temperature of 
HLX-1 in the soft state roughly follow the L oc corre¬ 
lation, consistent with the prediction of the standard disk 
model and with what is observed in stellar-mass XRBs 
(Servillat et al. 2011). In addition, the larger radius and 
lower temperature (compared with stellar-mass XRBs) 
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Figure 11. Long-term X-ray light curve of H1743—322 measured by Swift/HAT. It shows repeated outbursts nearly each year after MJD 
54700, similar to the outburst sequence of HLX-1. 
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Figure 12. Evolutions of the waiting time and recurrence time 
for the most recent outbursts of H1743—322. The zero outburst 
corresponds to the outburst peaking on about MJD 54470. The 
^Waiting, 'Flrecurrence and T2recurrence are defined in the Same way 
to HLX-1. 


are self-consistent with each other and with the increase 
in luminosity: that is, both scale as expected for a stan¬ 
dard disk when the BH mass is increased by a factor of 
^1000 to explain the correspondingly higher luminosity 
of the soft state. We conclude that all spectral evidence is 


consistent with the standard disk scenario, that is, with 
sub-Eddington, not strongly beamed emission, in which 
case HLX-1 must be an IMBH. 

Disk emission is not the only way to produce a ther¬ 
mal X-ray spectrum. Optically thick outflows caused 
by super-Eddington accretion may also produce a quasi- 
thermal X-ray spectrum with a low temperature and 
large radius, and high apparent (non-isotropic) luminos¬ 
ity if observed face-on (King & Pounds 2003; Begehnan 
et al. 2006; Fabrika et al. 2007; Poutanen et al. 2007; 
King 2009). The reason for this is that an optically thick 
outflow downscatters photons and lets them escape only 
through a collimating funnel along the polar axis (geo¬ 
metric beaming). The beaming factor b in the apparent 
luminosity (L = bLi^t, where Lint is the intrinsic lumi¬ 
nosity) for a face-on observer depends on the opening 
angle 6 of the funnel: for a small angle, b « 2/0^. Fol¬ 
lowing the model of King (2009), b ^ j73 and there¬ 
fore 9 ~ 12/m rad, where m = O.lc^M/LEdd- The in¬ 
trinsic luminosity for a super-Eddington accretor scales 
as Lint ~ LEdd(l + 0.6Inm) (Shakura & Sunyaev 1973; 
Poutanen et al. 2007). Combining the two factors, a 
face-on observer would measure an apparent luminosity 
L ~ (LEdd/73) rh^ (1 -|- 0.61nm). Thus, if the beaming 
scale holds, a stellar-mass BH with M < 20Mq accret¬ 
ing at m ~ 100 would be seen at L ~ 10"^^ erg s“^ by 
an observer looking down the collimating funnel (King 
& Lasota 2014; Lasota et al. 2015). There is no phys¬ 
ical reason preventing m from reaching values ~100 or 
even much higher (Wiktorowicz et al. 2015). The ques¬ 
tion is whether the opening angle of the funnel continues 
to scale as m“^ at such high accretion rates. In the 
model of King (2009), m ~ 100 corresponds to a fun¬ 
nel opening angle 0 ~ 7°. This is narrower than the 
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opening angle predicted by MHD simulations for com¬ 
parably high accretion rates: for example, it is « 20° 
in Ohsuga & Mineshige (2011) and « 30° in Jiang et al. 
(2014). Those simulations suggest that at the highest ac¬ 
cretion rates, the opening angle and the geometric beam¬ 
ing factor have a weaker dependence on m or reach an 
asymptotic value. In addition, when Compton down¬ 
scattering is also properly accounted for, the luminosity 
of the photons emerging from the collimating funnel is 
less than expected simply from the beaming factor (that 
is, a larger part of the intrinsic luminosity is transferred 
to cooler electrons in the outflow). As a result, the mod¬ 
els of Ohsuga & Mineshige (2011) predict that the geo¬ 
metrically beamed apparent luminosity may reach only 
« 22LEdd; Kawashima et al. (2012) predict peak ap¬ 
parent luminosities « 3 x 10^° erg s“^ from a 10 Mq 
BH. Both values are much less than required to explain 
the high state of HLX-1. Therefore, we argue that the 
possibility of reaching apparent luminosities « 10^^ erg 
s“ ^ via geometric beaming of a super-Eddington stellar- 
mass BH is not ruled out but appears unlikely based on 
the MHD simulations {e.g. Ohsuga & Mineshige 2011; 
Kawashima et al. 2012; Jiang et al. 2014), and has never 
been observationally proven in any real system, either. 

A second difficulty of the hyper-accretion scenario is 
that the spectra of the photon emission predicted to 
emerge from the funnel (Kawashima et al. 2012), at least 
for accretion rates up to to = 100, do not look like the 0.2 
keV thermal emission seen in HLX-1. It is true that the 
emerging spectrum is predicted to soften at higher accre¬ 
tion rates, but it retains a dominant Comptonized com¬ 
ponent up to ~ 10 keV, with a roll-over at energies of > 5 
keV (Figures 2 and 4 in Kawashima et al. 2012), similar 
to the spectrum seen in “normal” ULXs (making them 
perfectly consistent with the super-Eddington scenario). 
To our knowledge, so far, simulations have not explored 
what happens to the spectrum at to 3> 100. Moreover, 
neither analytic slim-disk models {e.g. Watarai et al. 
2001) nor MHD simulations {e.g. Kawashima et al. 2012) 
predict a spectral state transition between a hard, power- 
law-like, moderately super-Eddington state, and a purely 
thermal, extremely super-Eddington regime, as would 
be the case in HLX-1. (Recall that if we interpret the 
high/soft state of HLX-1 as due to geometric beaming at 
TO « 100, then the low/hard state would also be super- 
Eddington, with TO « 30.) Again, there is no physical ar¬ 
gument ruling out this possibility, but it has never been 
observed in real systems or predicted in simulations. 

Another clue comes from the brightness of the opti¬ 
cal counterpart, in particular, for its near-UV and blue 
emission. Early suggestions that this emission was from a 
massive, young star cluster (Farrell et al. 2012) have been 
disproved after it was found that the optical emission 
varies substantially between high and low X-ray states 
(Webb et al. 2014). For example, in the near-UV and U 
bands, the optical flux is « 1.5 mag higher in outburst 
than in quiescence (R. Soria et al. 2015, in prepara¬ 
tion). Thus, the most plausible interpretation is that 
at least the near-UV and blue emission are dominated 
by an irradiated accretion disk (Soria et al. 2012). The 
need for irradiation is crucial: a simple, non-irradiated 
Shakura-Sunyaev disk is too cold and dim to fit the ob¬ 
served optical/UV spectrum. The disk must intercept 


and reprocess an X-ray luminosity of a few times 10^® 
erg s“^ (Soria et al. 2012; Farrell et al. 2014). For any 
plausible solid angle subtended by the disk, this requires 
that the intrinsic X-ray luminosity in outburst be «10^^ 
erg s“^, and that it be approximately isotropic. Very 
few irradiation effects on the disk would be expected if 
the observed X-ray flux emerged from a collimated fun¬ 
nel perpendicular to the disk plane and was confined by 
an optically thick outflow. 

Taking into account all of the arguments and obser¬ 
vational constraints outlined above, we conclude that 
at this stage, the sub-Eddington IMBH scenario with 
canonical state transitions similar to the Galactic XRBs 
is the least implausible model. 

3.3. Disk Instability Model and Comparison with 
Galactic LMXBTs 

One of the objectives of this work is to quantify to what 
degree HLX-1 outbursts are similar to Galactic LMXBTs 
outbursts. The outbursts of bright Galactic LMXBTs 
are usually modeled as thermal-viscous disk instabilities 
(Lasota 2001, for a review). The correlation between 
the total radiated energy and the e-folding rise or decay 
timescale, and between the total radiated energy and the 
peak luminosity, found in Galactic LMXBTs (Yan & Yu 
2015), are consistent with the predictions of this model. 
Instead, the e-folding rise and decay timescales of HLX- 
1 are about two orders of magnitude smaller than the 
extrapolation of those correlations to a distance of 95 
Mpc (Figure 10). This is consistent with the argument 
(Lasota et al. 2011, 2015) that the thermal-viscous disk 
instability model cannot explain the short outburst du¬ 
ration in HLX-1. 

For the sake of argument, let us assume instead that 
we do not have any information on the distance to HLX- 
1, and let us try to constrain its distance according to 
the best-fitting results (log E = A + B x log r) of the 
E-t correlations in bright Galactic LMXBTs (Yan & Yu 
2015). We would obtain a distance of 540 kpc. Since 
the intrinsic scatters of those E-t correlations are quite 
large (r:! 0.4 dex), the uncertainty of the inferred distance 
can be up to a factor of 1.6. Other correlations we found 
in Galactic LMXBTs would be satisfied if HLX-1 were at 
a distance of 540 kpc (Figure 8, Figure 9 and Figure 10). 
With such a distance, we estimated the orbital period 
of HLX-1 («20 hr) according to the correlation between 
optical luminosity and X-ray luminosity and the orbital 
period found in Galactic LMXBs (van Paradijs & Mc- 
Clintock 1994), where the V-band magnitude of HLX-1 
near the 2012 outburst peak is from Webb et al. (2014), 
and we estimated the mass transfer rate according to 
the method in Coriat et al. (2012). The mass transfer 
rate and orbital period satisfy the criteria for a transient 
LMXB predicted by the thermal-viscous disk instability 
model (van Paradijs 1996; Coriat et al. 2012). In other 
words, the temporal properties of the HLX-1 outbursts 
are consistent with those of a stellar-mass X-ray tran¬ 
sient at several hundreds of kiloparsecs, which can be in¬ 
terpreted by the thermal-viscous disk instability model. 
Lasota et al. (2015) derived a relation between the peak 
luminosity and the decay timescale from the thermal- 
viscous disk instability model. According to the relation, 
the outbursts of HLX-1 with IMBH at 95 Mpc cannot be 
triggered by the thermal-viscous disk instability, and the 


11 


thermal-viscous disk instability model requires a stellar- 
mass accretor at a distance of less than 1 Mpc. 

This does not mean that HLX-1 is indistinguishable 
from a stellar-mass LMXBTs hypothetically located at 
that distance: other outburst properties of HLX-1 are 
very different. For example, the inner-disk tempera¬ 
ture of HLX-1 in the soft state is ftiO.2 keV (a distance- 
independent measurement), empirically and theoreti¬ 
cally inconsistent with the disk temperatures of Galactic 
LMXBTs in the same state {kT ~ 1 keV). Moreover, the 
fitted inner-disk radius of HLX-1 at 95 Mpc is R:!50,000- 
100,000 km depending on inclination angle, BH spin, and 
other detailed disk model assumptions. If HLX-1 were 
at a distance of 540 kpc, then the fitted inner-disk ra¬ 
dius would be R:i250-500 km. Once again, this value is 
inconsistent with the typical inner-disk radii of Galac¬ 
tic NS (i?i„ « 15 km) and BH LMXBTs (i?i„ « 50-100 
km). It is still consistent with the inner-most stable or¬ 
bit of more massive stellar BHs (M « 30-60Mq). How¬ 
ever, if the BH were in such mass range, and HLX-1 was 
at ~540 kpc, then its bolometric transition luminosity 
Lsoit-to-hardwouldbe ~ 10“^LEdd- This is inconsistent 
with the canonical state transition luminosity of BHs, 

for which Lsoft-to-hard ~ 1% - 4%LEdd- 

In summary, although the outburst timescales are sur¬ 
prisingly short for an IMBH at 95 Mpc (Lasota et al. 
2011, 2015), and would instead indicate a stellar-mass 
source at 540 kpc, this is not a viable solution for HLX-1 
because at that lower distance, its luminosity and peak 
temperature would no longer be self-consistent for the 
disk-dominated high/soft state. In addition, we do have 
other independent evidence that the distance to HLX-1 
is «100 Mpc: an optical emission line consistent with 
Ha redshifted to ~ 7000 km s“^, and its projected lo¬ 
cation in the sky in the middle of a cluster of galaxies 
at the same distance (Wiersema et al. 2010; Soria et al. 
2013). Gonversely, we have indirect evidence to disfa¬ 
vor a location at 540 kpc, because there are no known 
Local Group stellar structures (dwarf or satellite galax¬ 
ies, or star clusters) in that direction of the sky at that 
distance from the Milky Way. Therefore, we conclude 
that the anomalously short outburst timescales indicate 
that either the outbursts of HLX-1 with IMBH are not 
triggered by the thermal-viscous disk instability, which 
is consistent with the conclusions in Lasota et al. (2011, 
2015), or that the efficiency of the angular momentum 
transport in such a disk would have to be two orders of 
magnitude greater than in a standard disk, or that the 
circularization radius and outer radius of the accretion 
disk are two orders of magnitude smaller than naively 
estimated based on the size of the primary Roche lobe in 
a circular binary system. 

3.4. Comparison with Galactic Transient Be/XRBs 

High mass donor stars have been confirmed in some 
ULXs such as NGG 7793 P13 (Motch et al. 2011, 2014) 
and MlOl ULX-1 (Liu 2009; Liu et al. 2013). Therefore, 
it is also interesting to compare the temporal and spectral 
outburst properties of HLX-1 with those of Galactic tran¬ 
sient high mass X-ray binaries (HMXBs). Be/XRBs are 
the most common type of transient HMXBs: they consist 
of a NS bound to a Be star in an eccentric orbit. This 
kind of systems usually shows two types of X-ray out¬ 
bursts (Paul & Naik 2011; Reig 2011, for reviews). Type 


I outbursts are quasi-periodic, and likely triggered by 
the periastron passage of the NS. The amplitude of this 
type of outburst is smaller and the duration is shorter 
than for the HLX-1 outbursts; besides, we have already 
noted that the periastron-passage model for HLX-1 has 
serious difficulties. Type H outbursts usually show much 
larger amplitudes and longer durations, comparable with 
the outbursts of HLX-1. It has been suggested (Reig 
et al. 2006; Reig 2008; Becker et al. 2012; Reig & Nespoli 
2013) that such outbursts include a transition between 
distinct accretion states. In an HID, transient HMXBs 
move between different states either along a horizontal 
branch, at low luminosities, or along a diagonal branch, 
at higher luminosities (Reig & Nespoli 2013). On the 
diagonal branch, softer colors correspond to higher lu¬ 
minosities, as in Galactic LMXBTs and HLX-1. The 
physical mechanism is still unclear, but it is not due to 
a disk instability; more likely, it depends on how the ac¬ 
cretion flow is decelerated in the accretion column, and 
therefore on the balance between the accretion rate and 
magnetic field of the accreting NS (Reig & Nespoli 2013). 
As such, any phenomenological similarities between such 
outbursts and those of HLX-1 are most likely a coinci¬ 
dence. 

4. GONCLUSIONS 

We studied the X-ray outburst properties of HLX- 
1, using the full set of Swift/XRT monitoring observa¬ 
tions available to-date, which cover six outbursts from 
2009 to 2015 (plus a few sparse data points for the 
2008 outburst). HIDs confirm the presence of spectral 
state transitions during the outbursts (Godet et al. 2009; 
Servillat et al. 2011), similar to the cycle of low/hard 
and high/soft states in Galactic XRBs. We found 
that HLX-1 follows the nearly linear correlation between 
Thard-to-soft and Apeak (Figure 5) and the nearly lin¬ 
ear correlation between Lj-ise and Apeak (upper panel of 
Figure 8), and that the rise timescale is similar to those 
observed in stellar-mass Galactic sources; as a result, the 
hard-to-soft state transition also occurs at similar times 
after the beginning of an outburst. The similarity of the 
rise and decay timescales and outburst duration between 
HLX-1 at 95 Mpc and Galactic stellar-mass BHs remains 
an unexplained puzzle: on the one hand, it suggests that 
the characteristic size of the accretion flow is similar; on 
the other hand, we find no reason to dispute that the ac¬ 
cretion rate and the BH mass (and hence, the size of the 
ISGO) of HLX-1 are three orders of magnitude higher 
than in Galactic BHs. 

For Galactic LMXBTs, the correlations between 
the radiated energy and the e-folding rise and decay 
timescales, and between the radiated energy and the 
peak luminosity (Yan & Yu 2015) quantitatively sup¬ 
port the thermal-viscous disk instability model. HLX-1 
does not follow those correlations, which suggests that its 
outbursts are not triggered by the thermal-viscous disk 
instability (Lasota et al. 2011, 2015). If they were, the 
outer radius of the accretion disk would have to be two 
orders of magnitude larger than previously estimated, 
in order for a cold neutral zone to be present (Lasota 
et al. 2011); but in that case, the viscosity parameter 
would have to be two orders of magnitude higher than in 
Galactic LMXBTs, to obtain the same (short) rise and 


12 


decay timescales, and yet such high viscosity seems not 
plausible in the reality (King et al. 2007). 

We quantified the evolutionary trends of the main out¬ 
burst parameters of HLX-1 between 2009 and 2015. The 
outburst duration, total radiated energy, and e-folding 
decay timescale decrease along the outburst sequence, 
confirming previous results (Godet et al. 2013; Miller 
et al. 2014). The waiting time, recurrence time, and 
e-folding rise timescale increase with the outburst se¬ 
quence. The increase in the recurrence time rules out 
earlier suggestions that the recurrence time corresponded 
to the binary period, or indicates that the binary period 
increases in some extreme cases {e.g. Godet et al. 2014, 
2015). The increase of the waiting time and the decrease 
of the outburst duration show that HLX-1 spends more 
and more time in the quiescence phase and less and less 
time in the outburst phase; however, its peak luminosity 
continues to reach similar values of « 10"^^ erg s“^. 

This work was supported in part by the National 
Natural Science Foundation of Ghina under grant No. 
11333005 and 11350110498, by Strategic Priority Re¬ 
search Program “The Emergence of Gosmological Struc¬ 
tures” under grant No. XDB09000000, and the XTP 
project under grant No. XDA04060604, by the Shanghai 
Astronomical Observatory Key Project. Z.Y. acknowl¬ 
edges the support from the Knowledge Innovation Pro¬ 
gram of the Chinese Academy of Sciences and National 
Natural Science Foundation of China under grant No. 
11403074. R.S. acknowledges an Australian Research 
Councils Discovery Projects funding scheme (project No. 
DP 120102393). D.A. acknowledges support from the 
Royal Society. This research made use of Astropy, a 
community-developed core Python package for Astron¬ 
omy (Astropy Collaboration et al. 2013). This work made 
use of data supplied by the UK Swift Science Data Centre 
at the University of Leicester. This research has made 
use of the MAXI data provided by RIKEN, JAXA and 
the MAXI team, and the ASM data provided by the 
ASM/RXTE team. 

REFERENCES 

Arnaud, K. A. 1996, in Astronomical Society of the Pacific 
Conference Series, Vol. 101, Astronomical Data Analysis 
Software and Systems V, ed. G. H. Jacoby Sz J. Barnes, 17 [2.1] 
Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 
2013, A(^A, 558, A33 [4] 

Bachetti, M., Harrison, F. A., Walton, D. J., et al. 2014, Nature, 
514, 202 [1] 

Becker, P. A., Klochkov, D., Schonherr, G., et al. 2012, ASzA, 

544, A123 [3.4] 

Begelman, M. C., King, A. R., Sz Pringle, J. E. 2006, MNRAS, 
370, 399 [3.2] 

Belloni, T., Parolin, L, Del Santo, M., et al. 2006, MNRAS, 367, 
1113 [2.1] 

Belloni, T. M. 2010, in Lecture Notes in Physics, Berlin Springer 
Verlag, Vol. 794, Lecture Notes in Physics, Berlin Springer 
Verlag, ed. T. Belloni, 53 [2.1] 

Cash, W. 1979, ApJ, 228, 939 [2.1] 

Chen, W., Shrader, C. R., Livio, M. 1997, ApJ, 491, 312 [1, 
2.3] 

Coriat, M., Fender, R. P., Sz Dubus, G. 2012, MNRAS, 424, 1991 
[3.3] 

Davis, S. W., Narayan, R., Zhu, Y., et al. 2011, ApJ, 734, 111 
[2.1, 3.1] 

Dunn, R. J. H., Fender, R. P., Kording, E. G., Belloni, T., Sz 
Gabanac, C. 2010, Mon. Not. R. Astron. Soc., 403, 61 [2.1] 
Evans, P. A., Beardmore, A. P., Page, K. L., et al. 2007, A&A, 
469, 379 [2.1, 2] 

—. 2009, MNRAS, 397, 1177 [2.1, 2, 2.1] 


Fabrika, S. N., Abolmasov, P. K., &; Karpov, S. 2007, in lAU 
Symposium, Vol. 238, lAU Symposium, ed. V. Karas & 

G. Matt, 225-228 [3.2] 

Farrell, S. A., Servillat, M., Oates, S. R., et al. 2010, X-ray 
Astronomy 2009; Present Status, Multi-Wavelength Approach 
and Future Perspectives, 1248, 93 [3.y 
Farrell, S. A., Webb, N. A., Barret, D., Godet, O., Sz Rodrigues, 

J. M. 2009, Nature, 460, 73 [1, 2.1] 

Farrell, S. A., Servillat, M., Pforr, J., et al. 2012, ApJ, 747, L13 
[3.2] 

Farrell, S. A., Servillat, M., Gladstone, J. G., et al. 2014, 

MNRAS, 437, 1208 [3.2] 

Fender, R. P., Belloni, T. M., Sz Gallo, E. 2004, MNRAS, 355, 
1105 [1] 

Feng, H., Sc Soria, R. 2011, New A Rev., 55, 166 [1] 

Godet, O., Barret, D., Webb, N. A., Farrell, S. A., Sc Gehrels, N. 

2009, ApJ, 705, L109 [1, 2.1, 2.1, 3.2, 4] 

Godet, O., Lombardi, J. C., Antonini, F., et al. 2014, ApJ, 793, 
105 [1, 1, 4, 7, 4] 

Godet, O., Webb, N., Barret, D., et al. 2013, The Astronomer’s 
Telegram, 5439, 1 [1, 1, 4, 4] 

Godet, O., Webb, N., Barret, D., Sc Lombardi, J. 2015, The 
Astronomer’s Telegram, 6929, 1 [1, 7, 4] 

Godet, O., Plazolles, B., Kawaguchi, T., et al. 2012, ApJ, 752, 34 

[ 3-^1 . . 

Homan, J., Wijnands, R., van der Klis, M., et al. 2001, ApJS, 

132, 377 [2.1] 

Jiang, Y.-F., Stone, J. M., Sz Davis, S. W. 2014, ApJ, 796, 106 

[ 3-21 

Kalemci, E., Dinger, T., Tomsick, J. A., et al. 2013, ApJ, 779, 95 
[3.1] 

Kawashima, T., Ohsuga, K., Mineshige, S., et al. 2012, ApJ, 752, 
18 [3.2] 

Kelly, B. C. 2007, ApJ, 665, 1489 [5, 2.4] 

King, A., & Lasota, J.-P. 2014, MNRAS, 444, L30 [1, 3.2] 

King, A. R. 2009, MNRAS, 393, L41 [3.2] 

King, A. R., &; Pounds, K. A. 2003, MNRAS, 345, 657 [3.2] 

King, A. R., Pringle, J. E., & Livio, M. 2007, MNRAS, 376, 1740 
[4] 

Kong, A. K. H. 2011, ArXiv e-prints, arXiv:1102.3212 [1] 

Kong, A. K. H., Soria, R., & Farrell, S. 2015, The Astronomer’s 
Telegram, 6916, 1 [1] 

Kording, E., Rupen, M., Knigge, C., et al. 2008, Science (80-. )., 
320, 1318 [2.1] 

Krimm, H. A., Holland, S. T., Corbet, R. H. D., et al. 2013, 

ApJS, 209, 14 [2.5] 

Lasota, J.-P. 2001, New A Rev., 45, 449 [1, 3.3] 

Lasota, J.-P., Alexander, T., Dubus, G., et al. 2011, ApJ, 735, 89 
[1, 3.3, 4] 

Lasota, J.-P., King, A. R., & Dubus, G. 2015, ApJ, 801, L4 [1, 
2.4, 3.2, 3.3, 4] 

Liu, J. 2009, ApJ, 704, 1628 [3.4] 

Liu, J.-F., Bregman, J. N., Bai, Y., Justham, S., & Crowther, P. 

2013, Nature, 503, 500 [1, 2.1, 3.4] 

Maccarone, T. J. 2003, A&A, 409, 697 [3.1] 

Maccarone, T. J., & Coppi, P. S. 2003, Mon. Not. R. Astron. 

Soc., 338, 189 [2.1] 

Matsuoka, M., Kawasaki, K., Ueno, S., et al. 2009, PASJ, 61, 999 
[2.5] 

Miller, M. C., Farrell, S. A., Sz Maccarone, T. J. 2014, ApJ, 788, 
116 [1, 3, 5, 4] 

Miyamoto, S., Kitamoto, S., Hayashida, K., Sc Egoshi, W. 1995, 
ApJ, 442, L13 [2.1] 

Motch, C., Pakull, M. W., Grise, F., Sc Soria, R. 2011, 
Astronomische Nachrichten, 332, 367 [3.4] 

Motch, C., Pakull, M. W., Soria, R., Grise, F., Sz Pietrzyhski, G. 

2014, Nature, 514, 198 [1, 3.4] 

Munoz-Darias, T., Fender, R. P., Motta, S. E., Sz Belloni, T. M. 

2014, MNRAS, 443, 3270 [2.1] 

Ohsuga, K., Sz Mineshige, S. 2011, ApJ, 736, 2 [3.2] 

Paul, B., Sz Naik, S. 2011, Bulletin of the Astronomical Society of 
India, 39, 429 [3.4] 

Poutanen, j., Lipunova, G., Fabrika, S., Butkevich, A. G., Sz 
Abolmasov, P. 2007, MNRAS, 377, 1187 [2.1, 3.2] 

Reig, P. 2008, A&A, 489, 725 [3.4] 

—. 2011, Ap&SS, 332, 1 [3.4] 

Reig, P., Martmez-Nunez, S., Sz Reglero, V. 2006, A&A, 449, 703 
[3.4] 

Reig, P., & Nespoli, E. 2013, A&A, 551, Al [3.4] 

Servillat, M., Farrell, S. A., Lin, D., et al. 2011, ApJ, 743, 6 [1, 
2.1, 2.1, 3.1, 3.2, 4] 

Shakura, N. L, & Sunyaev, R. A. 1973, Astron. Astrophys., 24, 
337 [3.2] 


13 


Shen, R.-F., Barniol Duran, R., Nakar, E., & Piran, T. 2015, 
MNRAS, 447, L60 [2.1] 

Soria, R. 2013, MNRAS, 428, 1944 [2.3] 

Soria, R., Hakala, P. J., Hau, G. K. T., Gladstone, J. G., Sz Kong, 
A. K. H. 2012, MNRAS, 420, 3599 [3.2] 

Soria, R., Hau, G. K. T., & Pakull, M. W. 2013, ApJ, 768, L22 
[1, 3.3] 

Tang, J., Yu, W.-F., Yan, Z. 2011, Res. Astron. Astrophys., 11, 
434 [2.2] 

van Paradijs, J. 1996, Astrophys. J. Lett., 464, L139 [3.3] 
van Paradijs, J., &; McClintock, J. E. 1994, A&A, 290, 133 [3.3] 
Watarai, K.-y., Mizuno, T., S>z, Mineshige, S. 2001, ApJ, 549, L77 
[3.2] 


Webb, N., Cseh, D., Lenc, E., et al. 2012, Science, 337, 554 [1, 
2.1, 3.1] 

Webb, N. A., Barret, D., Godet, O., et al. 2010, ApJ, 712, L107 
[ 2-11 . 

Webb, N. A., Godet, O., Wiersema, K., et al. 2014, ApJ, 780, L9 
[3.2, 3.3] 

Wiersema, K., Farrell, S. A., Webb, N. A., et al. 2010, ApJ, 721, 
L102 [1, 3.^ 

Wiktorowicz, G., Sobolewska, M., Sadowski, A., Sz Belczynski, K. 

2015, ArXiv e-prints, arXiv:1503.08745 [3.2] 

Yan, Z., Sz Yu, W. 2010, Sci. China G Phys. Astron., 53, 147 [2.1] 
Yan, Z., & Yu, W. 2015, ApJ, 805, 87 [1, 2.3, 2.4, 2.4, 3.3, 4] 

Yu, W., & Dolence, J. 2007, ApJ, 667, 1043 [2.2] 

Yu, W., van der Klis, M., Sz Fender, R. 2004, ApJ, 611, L121 [2.2] 
Yu, W., Sz Yan, Z. 2009, ApJ, 701, 1940 [2.1, 2.2] 


